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ABSTRACT 

General criteria to check the positivity of the distribution function (phase-space consistency) 
of stellar systems of assigned density and anisotropy profile are useful starting points in Jeans- 
based modeling. Here we substantially extend previous results, and we present the inversion 
formula and the analytical necessary and sufficient conditions for phase-space consistency 
of the family of multi-component Cuddeford spherical systems: the distribution function of 
each density component of these systems is defined as the sum of an arbitrary number of 
Cuddeford distribution functions with arbitrary values of the anisotropy radius, but identical 
angular momentum exponent. The radial trend of anisotropy that can be realized by these 
models is therefore very general. As a surprising by-product of our study, we found that the 
"central cusp-anisotropy theorem" (a necessary condition for consistency relating the values 
of the central density slope and of the anisotropy parameter) holds not only at the center, 
but at all radii in consistent multi-component generalized Cuddeford systems. This last result 
suggests that the so-called mass-anisotropy degeneracy could be less severe than what is 
sometimes feared. 
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1 INTRODUCTION 

In the study of stellar systems based on the "p-to-/" approach 
(where p is the material density and / is the associated phase- 
space distribution function, hereafter DF; e.g. see Bertin 2000, 
Binney & Tremaine 2008), the density distribution is given, and 
specific assumptions on the internal dynamics of the model are 
made. In some special cases inversion formulae exist so that the 
DF can be obtained, usually in integral form or as series expan- 
sion (see, e.g., Fricke 1952; Lynden-Bell 1962; Osipkov 1979; 
Merritt 1985; Dejonghe 1986, 1987; Cuddeford 1991; Hunter & 
Qian 1993; Ciotti & Bertin 2005). Once the DF of the system is 
derived, a non-negativity check is (or should be) performed, and 
in case of negative values the model must be discarded as unphys- 
ical. Indeed, a minimal but essential requirement to be met by the 
DF (of each component) of a stellar dynamical model is positivity 
over the accessible phase-space. This requirement, the so-called 
phase-space consistency, is much weaker than the model stability, 
but it is stronger than the fact that the Jeans equations have a phys- 
ically acceptable solution. However, the difficulties inherent in the 
operation of recovering analytically the DF prevent in general a 
simple consistency analysis, and numerical inspection of the inver- 
sion integral is required. As a consequence, the reasons underlying 
consistency or inconsistency of a proposed model are somewhat 
obscured by the numerical nature of the solution. Fortunately, cri- 
teria for phase-space consistency that can be applied without an 



explicit recovery of the DF are known and widely used. For exam- 
ple, analytical necessary and sufficient conditions for consistency 
of multi-component systems with Osipkov-Merritt anisotropy (Os- 
ipkov 1979, Merritt 1985, hereafter OM) were derived in Ciotti & 
Pellegrini (1992, hereafter CP92; see also Tremaine et al. 1994) and 
applied in several investigations (e.g., Ciotti 1996, hereafter C96; 
Ciotti 1999, hereafter C99; Ciotti & Lanzoni 1997; Ciotti & Mor- 
ganti 2009, hereafter CM09; Ciotti, Morganti & de Zeeuw 2009). 
Such conditions revealed not only simple and useful to investigate 
the phase-space consistency of OM models, but also helpful to elu- 
cidate the different roles of total potential, orbital anisotropy, and 
stellar and dark matter density profiles in making a model unphys- 
ical. 

More recently, the "central cusp-anisotropy theorem" (An & 
Evans 2006, hereafter AE06), a necessary condition for consistency 
relating the values of the central density slope and of the anisotropy 
parameter (3 (see equation (4)) has been provecQ. This condition 
was derived for constant anisotropy systems, and then generalized 
asymptotically to the central regions of spherical systems with arbi- 
trary anisotropy distribution. A remarkable property of the density 
slope-anisotropy inequality is that it actually holds rigorously at ev- 
ery radius in constant anisotropy systems, and not only at their cen- 
ter (AE06). Surprisingly, in CM09 we showed that the CP92 neces- 
sary condition for model consistency can be formally rewritten as 
the AE06 inequality, that consequently holds at each radius not only 
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1 The same inequality was also reported in equation (28) in de Bruijne et 
al. (1996). 
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in constant anisotropy systems, but also in multi-component OM 
systems! This curious result prompted us to investigate the phase- 
space consistency of Cuddeford (1991) anisotropic systems, as they 
generalize both constant and OM anisotropy and an explicit inver- 
sion formula exists, so that necessary and sufficient conditions for 
consistency can hopefully be found, extending those of CP92. In 
addition, Cuddeford anisotropy allows to explore systems in which 
the central regions may be tangentially anisotropic, at variance with 
the OM cases. 

Actually, we found it possible to extend our study to the very 
general case of multi-component, generalized Cuddeford systems, 
i.e. spherical systems in which the DF of each distinct density com- 
ponent is assumed to be the sum of an arbitrary number of Cud- 
deford DFs with arbitrarily different anisotropy radii, but identi- 
cal angular momentum exponent (see equation 1221 ). In this pa- 
per we show how the family of necessary and sufficient conditions 
for model consistency can be derived for generalized Cuddeford 
anisotropic systems. We also found that the first of the necessary 
conditions coincides again with the density slope-anisotropy the- 
orem, thus demonstrating that such inequality must be satisfied at 
all radii also in the whole family of consistent, multi-component 
generalized Cuddeford systems. 

The paper is organized as follows. In Section 2 we recall 
the fundamental properties of OM and constant anisotropy sys- 
tems, and the associated consistency criteria. In Section 3 we de- 
rive the family of consistency criteria for the larger class of multi- 
component galaxy models with generalized Cuddeford anisotropy. 
Then, in Section 4 some illustrative applications of the new phase- 
space consistency criteria are presented, and in Section 5 the main 
results are summarized, with a brief discussion of the relevance of 
the new findings for the mass-anisotropy degeneracy problem. In 
the Appendix we prove that the first of the necessary conditions for 
phase-space consistency of multi-component generalized Cudde- 
ford systems can be rewritten as the density slope-anisotropy in- 
equality, that must hold at all radii. 



2 CONSISTENCY CRITERIA FOR 

MULTI-COMPONENT OSIPKOV-MERRITT 
SYSTEMS 

In this Section we summarize the main features of the OM inversion 
procedure, focusing on the arguments upon which the derivation 
of the CP92 necessary and sufficient conditions for phase-space 
consistency is based: in fact, similar arguments will be applied to 
generalized Cuddeford systems in Section 3. 

To fix the notation, we say that a multi-component stellar sys- 
tem described by a sum of different density components pi is called 
consistent if each DF fi is non-negative over the whole accessible 
phase-space. However, as all the conditions presented in this pa- 
per hold for each pi, for simplicity from now on the index i is not 
indicated, except when required for clarity. 

The OM prescription assumes that the DF supporting each 
density component depends on the energy and on the angular mo- 
mentum modulus of stellar orbits as 



f = f(Q), 



j 2 - 



(1) 



and f(Q) = for Q < 0. In the formula above E = * T - 
v 2 /2 is the binding energy per unit mass, \1/t = — $t, where <E>t 
is the potential due to the combined effect of all the components 
Pi, and r a is the so-called anisotropy radius of each component 



(e.g. see Binney & Tremaine 2008). Each density component of 
a multi-component OM system is characterized by a DF of the 
family QJ, in general with different r a i: therefore, unless all the r a i 
are identical, a multi-component OM system is not an OM system. 
It is easy to prove that the DF of each component is related to its 
density profile as 

/d 3 « = 2%/87ryl(r,r a ) J V*t - Qf(Q)dQ, (2) 

where 



ri + r 2 



(3) 



The radial dependence of the associated anisotropy parameter, a 
quantity designed to measure the differences between the tangential 
(of ) and radial (of) velocity dispersions, is 



r 1 + ri 



(4) 



(Merritt 1985), so that the orbital distribution is isotropic at the 
center and increasingly radially anisotropic with radius. Note that 
consistency implies j3 < 1. With the introduction of the so-called 
"augmented density" 



g(r) 



A(r, r a ) 



p{r) 



(5) 



it is possible to recast equation (0 in a form suitable for Abel in- 
version, and after the differentiation one obtains 



dg 



= y/8n 



f(Q)dQ 
V*t -Q 



(6) 



where the function g is intended to be expressed in terms of ^x, by 
the elimination of radius. As first solved by Eddington (1916) for 
the isotropic case in which Q — £ , equation ^ can be inverted as 



f(Q) 



d 



dg 



d^q 



VStv 2 dQ J d* T VQ- *T 



d 2 g d^T 



(7) 



(Osipkov 1979), where the second identity above holds for untrun- 
cated systems with finite total mass. Equation (|6j is also of central 
importance in the derivation of the CP92 necessary condition: 

Theorem [CP92, CM09] A necessary condition (NC) for the 
non-negativity of the DF of each density component p in a multi- 
component OM system is 



dg 



> 0, 0<*<*(0), 



(8) 



where g is the augmented density in equation lO, and ^ is the rel- 
ative gravitational potential of the considered density component. 
The NC can be rewritten in terms of the logarithmic density slope 
7(r) = — d\np/d\nr and of the anisotropy parameter f3(r) as 



7 (r) > 2p(r), Vr. 



(9) 



In addition, a weak sufficient condition (WSC) for the non- 
negativity of each DF is 

d I dg 



( > 0, < * < *(0). 



(10) 



Proof: see CP92, C96, and CM09. Here we just recall that the NC 
is obtained by assuming a positive f(Q) in equation ((6), while 



Consistency of generalized Cuddeford systems 3 



the WSC by requiring the positivity of the integrand in the sec- 
ond equation 0, i.e. inequality dlOt is nothing else that a rewriting 
of d 2 g/d^ 2 T > 0. 

Of particular relevance for the following discussion is inequal- 
ity {9]l, an unexpected extension of the "central slope-anisotropy 
theorem": 

Theorem [AE06] In all consistent constant anisotropy sys- 
tems (with < 1 /2) necessarily 



7 (r) > 2(3, Vr. 



(11) 



Moreover, the same inequality holds asymptotically at the center 
(i.e., for r — + 0) of any consistent spherical system with generic 
anisotropy profile. 
Proof: see Section 2.1.1 in AE06. 

For completeness, we recall that systems with constant 
anisotropy are generated assuming a DF of the form 



/ = J 2a h(£ ), 



(12) 



where h(£) is a positive function, and a > —1 is a real number 
(see Section 3; see also Binney & Tremaine 2008). In such models 
the anisotropy parameter is 



(3(r) 



-a. 



(13) 



so that for a > they are characterized by tangential anisotropy, 
while for — 1 < a < the orbital anisotropy is radial. The proof 
of identity < | 1 3 1 > and the inversion formula analogous to are not 
given here, being obtained as special cases of the Cuddeford sys- 
tems described in the next Section. 



3 CONSISTENCY CRITERIA FOR 

MULTI-COMPONENT GENERALIZED CUDDEFORD 
SYSTEMS 

We begin this Section by recalling the main features of the inver- 
sion for Cuddeford (1991) systems. Then, in Section 3.2 the family 
of multi-component generalized Cuddeford systems is introduced 
and the inversion formula obtained, together with the associated 
consistency conditions. 



3.1 Cuddeford systems 

An interesting generalization of OM and constant anisotropy sys- 
tems was proposed by Cuddeford (1991; see also Ciotti 2000, 
Chapter 10) assuming 



/ = J 2a h(Q), 



(14) 



where a > —1 is a real number and Q is defined as in equation 0: 
isotropic models then correspond to a — and r a — > oo. Equa- 
tion d!4t can be used to describe both the OM models (for a = 0) 
and the constant anisotropy models (for r a — > oo). In particular, 
the anisotropy parameter takes now the simple form 



/8(r) = 



+ ri 



(15) 



(see equations (AT)-[A5]). Therefore, when a > the anisotropy is 
tangential in the inner regions where r < ^/ar^, and radial for r > 
y/ar a ; in the limit a — + oo, the orbital structure is fully tangentially 
anisotropic (i.e., (3 — > — oo). Instead, when — 1 < a < the 
models are radially anisotropic everywhere, independently of the 
value of r a ; moreover, in the limit a — * — 1 equation \15\ gives 
(3 — > 1, so that the velocity anisotropy is completely radial. 



The DF of a Cuddeford system and its spatial density are re- 
lated as 



p(r) = 2V8tt A(r,a) 
where 

A{r,a) = 2°" 



i^T -Q) a+1/2 h(Q)dQ, 



r(a + i) 



(16) 



(17) 



v T(a + 3/2) (l + r 2 /^) Q+1 ' 

and r(a;) = (x — 1)! is the gamma function. As expected, equa- 
tion ((2} is reobtained for a = 0, while the convergence of the an- 
gular part of the integral over the velocity space requires a > — 1. 
In analogy with the discussion in Section 2, the augmented density 

P 



A{r,a) 

2 1 "" r(q + 3/2) 
~W T(a + 1) 



1 + 



p(r) 



(18) 



is introduced, and a simple inversion formula, similar to equa- 
tion (0, permits to recover the DF from the density profile. In fact, 
after 



m — int 



hi) 



+ 1 



(19) 



differentiation^ with respect to \?t, equation J 16b can be Abel in- 
verted (Cuddeford 1991). In practice, one must perform enough 
differentiations as to produce a negative exponent (> —1) in the 
power-law kernel of integral | |16I >. 

When a > — 1 (i.e., m > 0) but a is not half-integer, 

<-l) m+1 cosmr r(a + 3/2 - m) 



h(Q) 



dQ 



T(a + 3/2) 



2y/Sn 2 

<! f Q cTg 

rf*™ (Q - * T ) a+3/2 - m 
cos an r(a + 3/2 — m) 



\m+l 



2VEtt 2 



T(a + 3/2) 



(Q " * T )«+3/2-m : 



(20) 



where the last identity holds for untruncated systems with finite 
total mass, and the OM inversion formula is reobtained for a = 
0. 

When a is half-integer, i.e. a = m — 3/2 and m = 1, 2, 
the solution of the Volterra equation d!6t is given by 



HQ) 



2V8Tv(r 



d m Q 



(21) 



*T = Q 



and the DF is recovered analytically avoiding integratior0. 



3.2 The consistency criteria and the density slope-anisotropy 
inequality 

As we now show, the inversion formulae l |20t and d21t still hold for 
the more general case of multi-component, generalized Cuddeford 
systems, in which the DF associated with each density component 

2 int(ic) means the largest integer < x. For example, int(l/2) = and so 
m = 1 for OM models. 

3 In equation (30) of Cuddeford (1991) the (m — 1)! at the denominator is 
missing. See also equations (49) and (51) of Baes & Dejonghe (2002). 
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is made by the sum of an arbitrary number of Cuddeford DFs with 
arbitrary positive weights Wi and possibly different anisotropy radii 
r a i (but same h function and angular momentum exponent), i.e. 



/ = J 2a ^ WjhjQj 



: £ 



2rl 



(22) 



The different density components of a multi-component general- 
ized Cuddeford system will have, in general, a different value of 
q and a different function h(Q). As should be clear, all the results 
presented in Sections 2 and 3.1 hold as special cases of the follow- 
ing treatment. 

Of course, the orbital anisotropy distribution characteristic of 
DF ( 122b is not a Cuddeford one: as shown in the Appendix, the 
anisotropy function /3(r) of each density component is given by 



,2 \a + 2 



0(r) 



(a + 1) 



£^/(l + r 2 /: 



(23) 



Quite general anisotropy profiles can be obtained by specific 
choices of the weights ra,, the anisotropy radii r a i, and the expo- 
nent a. However, near the center /3(r) ~ —a, and /3(r) ~ 1 for 
r — > oo, independently of the specific values of to, and r a i. 

We now show that an Abel inversion formula identical to equa- 
tions d20b-J21b can be found for a DF of the family J22b . In fact, it 
is immediate to verify that equation j 1 6b still holds, where now the 
radial function is 



A(r,a) 



r(q + i) 

2 r(a + 3/2) (l + r 2/ r 2.) Q+ i > 



\ " 



(24) 



and so, once the new augmented density g — p/A is defined, the 
function h in equation ( 122b can in principle be recovered. There- 
fore, it is obvious that the same arguments used to derive the neces- 
sary and sufficient conditions for consistency of OM models can be 
repeated also for each density component of multi-component gen- 
eralized Cuddeford systems. However, as m differentiations with 
respect to must be performed on the integral J 16b before the 
inversion, we now obtain m necessary conditions and a sufficient 
condition. Surprisingly, as in the case of OM models, we found that 
the first of the necessary conditions for consistency can be rewritten 
as the density slope-anisotropy theorem which must hold at every 
radius. These results are summarized in the following 

Theorem Each density component in a consistent multi- 
component generalized Cuddeford system with a not half-integer 
obeys m necessary conditions (NCk): 



*f >0, 



k = 1,2, 



(25) 



where m is given by equation l |19b . In particular, the NCi can be 
rewritten as the density slope-anisotropy inequality 



7 (r) > 2/3(r), Vr. 



(26) 



Moreover, a sufficient condition for the non-negativity of the DF 
of each component is 



> 0. 



(27) 



d m+1 g 

Proof: A proof of the m necessary conditions ( |25b is obtained 
by repeated differentiation of the augmented density g (see equa- 
tions 1 1 61 - 1181 , where now A(r, a) is given by equation 1241 ) with 
respect to the total potential ^x, and by the assumption that h(Q) 
is a positive function. The sufficient condition l !27b is derived just 



by imposing the positivity of the integrand in the second of identi- 
ties ( 120b . Finally, we refer to the Appendix for a proof of inequal- 
ity t26t . 

Of course, in the special cases of a = m — 3/2 and m = 
1, 2..., equation ( 121b provides, in addition to the m — 1 necessary 
conditions (25), the necessary and sufficient condition for consis- 
tency of the specific component, i.e. 



d" 



> 0. 



(28) 



As expected, the NC and the WSC derived in CP92 are reobtained 
as special cases of the new theorem for a = (i.e. m — 1). In 
applications, as those that will be presented in Section 4, it is useful 
to express equations ( 125b and ( 127b in terms of the radius. From the 
relation d^^/dr = -GMT(r)/r 2 , where Mt(t) is the total mass 
enclosed by the radius r, the NCi can be expressed as 



d /<o, 

dr 



(29) 



which is a second alternative formulation of the density slope- 
anisotropy theorem in addition to equation ( 126b . Following the 
same approach, the NC2 can be also expressed as 



d_ 

dr 



dp 



M T (r) dr 



> 0, 



(30) 



and so on, with the sign of the NCk inequality alternating with in- 
creasing k. Finally, note that NCi is the sole condition in which 
only the augmented density profile of the specific density compo- 
nent appears, while in the higher order NCk the total mass profile 
tfr(r) is also involved. 



4 SOME ILLUSTRATIVE CASES 

In the previous Section we derived the family of necessary and suf- 
ficient conditions for phase-space consistency of each density com- 
ponent of generalized Cuddeford systems, and we showed that the 
density slope-anisotropy inequality holds at every radius. 

We now present a simple application of the new consistency 
criteria, and we address two natural questions concerning phase- 
space consistency of Cuddeford systems. The first is related to the 
fact that for a > 1/2 we have, at variance with the OM case 
(a — 0), more than one necessary condition for the non-negativity 
of the DF. Which necessary condition is stronger? Or, more quan- 
titatively, which of the NCk gives a consistency limit closer to the 
true one (that would be derived from the DF)? The second ques- 
tion is: for a given density profile, what is the effect of tangential 
anisotropy on consistency? Will the minimum anisotropy radius 
increase or decrease at increasing a, i.e. at increasing tangential 
anisotropy? The set of necessary conditions and the dependence of 
their number on a through equation l |19b suggest a simple approach 
to address the two issues above. Consider an assigned density pro- 
file, representing a component in a multi-component generalized 
Cuddeford system: what is the behaviour of the consistency region 
in parameter space at increasing ctl At increasing a the number of 
necessary conditions increases: of course each additional necessary 
condition can only reduce the consistency region in the parame- 
ter space. In addition, when a increases so that m given by equa- 
tion l |19b increases by 1, the former sufficient condition NC m +i 
becomes the last of the necessary conditions for the new model. 

We now illustrate the procedure, and discuss the two questions 
presented above, by investigating the phase-space consistency of 
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Figure 1. Consistency limits on the normalized anisotropy radius s a = 
Ta./r c for the one-component 7 = model with Cuddeford anisotropy. 
The solid curves mark the limits imposed by NCi and NC2: models with 
the pair (a,s a ) in the shaded regions are certainly inconsistent; models 
above the dashed curve (where NC2 is the sufficient condition 1271 ) are 
certainly consistent. Solid dots are the true lower limits for s a derived from 
the DF. Recall that a = refers to the OM model. No consistent models 
exist for a > 3/2. 



the widely used 7-models (Dehnen 1993, Tremaine et al. 1994; see 
equation Bit below). We do this in the most simplified form, i.e. 
in the case of a one-component Cuddeford system; in other words, 
in equation d22b we restrict to i = 1. The detailed study of 7 = 0, 
7 = 1 (Hemquist 1990), and 7 = 2 (Jaffe 1983) models will also 
allow us to explore the combined effect of the inner density slope 
and of tangential anisotropy (Section 4.1), while the additional role 
played by the external density slope will be discussed in Section 4.2 
by using one-component n-j models with Cuddeford anisotropy 
(see equation J34t below). We recall that a consistency analysis of 
OM anisotropic 71-7 models was done in CM09. 



4.1 The one-component 7-models 

We start by considering the general 7-model, whose dimensionless 
density profile and mass enclosed inside radius r are given by 



p(r) 



ST(1 + S) 4 -T : 



M(r) = (^-) \ < 7 < 3, 



(3D 
(32) 



where s = r/r c is the radius normalized to the "core" radius r c . It 
is trivial to show that the consistency properties of one-component 
models are independent of the mass and density normalization 
scales. 

In the following we will study the NCk functions by using 
their radial formulation (equations 1291 - 1301 with the augmented 
density of equation 1181 ). Indeed, in common situations the elimi- 
nation of the radius from the density profile in favour of the gravi- 
tational potential, needed to evaluate equation \25\ . is not feasible. 
For this reason we prefer to study the consistency conditions by us- 



Figure 2. Consistency limits on the normalized anisotropy radius s a = 
Te>/r Q for the one-component Hernquist (7 = 1) model with Cuddeford 
anisotropy. Different curves have the same meaning as in Fig.Q] The dotted 
line connecting the solid dots has been obtained from Table 1 in Baes & 
Dejonghe (2002). 



0.3 




Figure 3. Consistency limits on the normalized anisotropy radius s a = 
r a /fc for the one-component Jaffe (7 = 2) model with Cuddeford 
anisotropy. Different curves have the same meaning as in Fig. [T] Note that 
the NCi coincides with the rr-axis. 



ing their radial expressions, as this procedure is always viable, for 
whatever density profile expressed as a function of radius. 

We begin by noticing that from AE06 theorem we already 
know that a > —7/2 is required at the center of the density 
distribution (31); this condition must be combined with a > — 1 
(see Section 3.1). Once the appropriate augmented density is de- 
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fined, the radial NCi for Cuddeford anisotropic 7-models read 
with equation l !18b reduces to 

sl[2s{2 + a) + 2a + 7 ] + s 2 (2s + 7 -2) > 0, Vs, (33) 

thus establishing a relation between a and the normalized 
anisotropy radius s a = r a /r c . 

As expected, for s — the inequality above reduces to the 
AE06 limitation. However, as the condition d33t must hold over the 
entire radial range, we can now derive limitations on the minimum 
allowed anisotropy radius s a as a function of 7 and a. The general 
formula is simple but here we prefer to focus on the special cases of 
7 = 0, 1, and 2. The NCi is represented by a solid curve in Figs[T] 
[2] [3] respectively for the 7 = 0, 1, and 2 models; of course, while 
a is restricted to positive values when considering the 7 = case, 
the q axis begins at a — —0.5 for 7 = 1 models, and finally the 
AE06 limitation in the 7 = 2 case is a > — 1 (coincident with the 
value required by convergence of the integral in equation 1161 ). 

In the three figures, all points below the solid NCi curve cor- 
respond to unphysical models, while points above may represent 
consistent models. The solid dots are the true lower limits on s a 
determined by direct inspection of the DF for representative val- 
ues of a: the 7 = 1 case was already given by Baes & Dejonghe 
(2002). Note that for the Jaffe model (Fig.0, the NCi actually co- 
incides with the abscissae axis, i.e. it is satisfied for all values of 
a > — 1, independently of the value of s a . In the case of 7 = 1 and 
7 = 2 models, when a — —1/2 equation d21t provides the DF, and 
so the true critical anisotropy radius s a can be easily determined 
(black dots in Fig.|2|and|3>. Of course, these values coincide with 
those obtained from the NCi for a = —1/2, as should be clear 
from the discussion in Section 3. The solid dots at a = represent 
instead the DF-derived lower limit for the minimum anisotropy ra- 
dius for the corresponding OM models (see, e.g., C96). As the NCi 
for —0.5 < a < 0.5 is just a necessary condition for consistency, 
while the NC2 provides a sufficient condition for consistency (i.e. 
all points above the NC2 dashed curve correspond to consistent 
models), it is not surprising that for all the three density models the 
DF-derived limit on the anisotropy radius in the OM case is con- 
tained in the region delimited by the NCi and the NC2 (see C99, 
Table 1). 

As we increase a, when we reach the value a — 1/2 the NC2 
function becomes the model DF, and so the DF-derived lower limit, 
represented with a black dot, coincides again with the critical curve. 

For 1/2 < a < 3/2, the NC2 becomes a new necessary con- 
dition, and therefore all points in the shaded area below the solid 
NC2 curves in the three figures correspond to unphysical models. 
Note how the NC2 provides more stringent limits than the NCi. 
Consistently with the nature of the NC2, the black dots represent- 
ing the limits on s a derived from the DF for a = 1 lie above the 
NC2 curve. Of course, in this range of values of a the NC3 is 
the sufficient condition for phase-space consistency. However, an 
asymptotic expansion of the NC3 for s — ■> 00 easily shows that this 
condition is violated, independently of the value of s a and a. This 
fact poses no problem in the range 1/2 < a < 3/2, as NC3 is a 
sufficient condition there, but as soon as a becomes larger than 3/2 
the NC3 becomes necessary, and the whole family of 7-models 
with Cuddeford anisotropy becomes inconsistent. We note that the 
q = 3/2 limitation was already determined by Baes & Dejonghe 
(2002) for Hernquist models with Cuddeford anisotropy. Quite sur- 
prisingly, by using the NC3 we found that the limitation a < 3/2 
holds for the entire family of 7-models, no matter which value of 
7 is considered. The reason is due to the fact that the external den- 
sity slope of 7-models is 4, independently of the value of 7. Thus, 



while the lower limit on a is due to the central density slope, the 
external density slope limits the amount of tangential anisotropy 
that can be supported by the models. This indication is very inter- 
esting, because it means that the external regions (where anisotropy 
is almost completely radial, see equation 1151 ) are able to affect the 
inner dynamics. We will discuss such issue in the next Section 4.2. 

As a final remark, we note that a comparison of Figs Q] [2] and 
[3] confirms qualitatively the trend already found in Carollo et al. 
(1995), C96 and C99 for one-component OM models. In practice, 
at fixed a the minimum anisotropy radius increases at decreasing 
inner density slope 7, i.e. centrally flatter density profiles are less 
able to sustain radial anisotropy than steeper density profiles, even 
in presence of a central tangential anisotropy. This is shown by the 
smaller and smaller shaded areas at fixed a and increasing 7, and 
by the corresponding smaller exact values of the minimum s a indi- 
cated by the solid dots. 



4.2 The effect of the external density slope 

As we have seen, no Cuddeford anisotropic 7-model exists for 
a > 3/2, and this independently of the value of 7. The fact that 
the critical upper limit of a is independent of 7 is a clear indication 
of the importance of the slope of the outer density profile on the 
central anisotropy. However, being the limit imposed on a, it im- 
plies that the central regions are the ones affected. A hint to under- 
stand this phenomenon is given by inspection of the three figures: 
in fact, note how the consistency region in the (a,s a ) space reduces 
at increasing a, in the sense that for increasing a the minimum 
value of s a increases. This means that when the central regions are 
forced to be more and more tangentially anisotropic, the external 
regions (where s a determines the amount of radial anisotropy) must 
be more and more isotropic. Therefore, we conclude that the origin 
of the inconsistency at high a is a combination of the forced tan- 
gential anisotropy and the radial orbits arriving from the external 
regions of the system. 

To better understand this behaviour, we now consider the one- 
component, Cuddeford anisotropic n-7 models, whose normalized 
density profile is given by 

p=— F —— , 0< 7 <3, n>3 (34) 

st(1 + s) n -T ~ 

(see CM09). If our previous argument is correct, then the maxi- 
mum value of a should increase at increasing n, as less and less 
mass is contained outside the core radius at increasing n, so that 
less and less radial orbits can affect the inner regions. Unfortu- 
nately, for generic (non integer) values of n the mass contained 
within r cannot be expressed in terms of elementary functions. 
However, it is possible to perform an asymptotic analysis at large 
radii of the NCk (with some care, as differentiation of asymp- 
totic relations is usually not permitted, e.g. see Bender & Orszag 
1978). The appropriate way to perform the analysis in this case is 
to use the NCk formulated in terms of the potential (equation 1251 ). 
and to adopt the asymptotic expansion for the relative potential 
\& = 1/s + 0(1/ s 2 ). Following this approach it can be proved 
that, independently of the value of the inner density slope 7, the 
critical value of a increases with n: for example, when n = 5 it is 
required that a < 5/2, when n — 6 that a < 7/2, and so on. This 
confirms the previous conjecture. 
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5 DISCUSSION AND CONCLUSIONS 

In a natural extension of previous investigations (CP92, C96, C99, 
AE06, CM09), we searched for phase-space consistency criteria 
for multi-component spherical systems. We found that inversion 
formulae and necessary and sufficient conditions for consistency 
can actually be derived for multi-component generalized Cudde- 
ford systems. Such systems contain as very special cases OM, con- 
stant anisotropy, and Cuddeford models. The main results of our 
study can be summarized as follows: 

(i) New phase-space consistency criteria, i.e. necessary and suf- 
ficient conditions for the DF non-negativity, are derived for multi- 
component, generalized Cuddeford systems. At variance with the 
simpler case of OM models, the presence of tangential anisotropy 
leads to a family of necessary conditions, that can be written as 
simple inequalities involving repeated differentiations of the aug- 
mented density expressed as a function of the total potential. 

(ii) It is shown that the first of the necessary conditions for con- 
sistency can be reformulated as the density slope-anisotropy the- 
orem, which therefore is proved to hold not only at the center but 
also at all radii for each density component of multi-component 
generalized Cuddeford models. 

(iii) The first necessary condition is the only condition indepen- 
dent of the other density components of the model. All the other 
(more stringent) conditions depend on the total density distribution 
of the model. 

(iv) All the conditions can be reformulated in term of the radius, 
so that they can be tested also for models in which the total potential 
cannot be expressed by using elementary functions, or when the 
radius cannot be eliminated in favour of the potential. 

(v) The new phase-space consistency criteria are applied to 
one-component 7-models with Cuddeford anisotropy. It is found 
that for increasing tangential anisotropy in the central regions the 
minimum anisotropy radius for consistency increases, i.e. the ex- 
ternal regions must be less and less radially anisotropic. No con- 
sistent 7-models exist for a > 3/2, independently of the value of 
the central density slope 7. Baes & Dejonghe (2002) already found 
this limitation by direct inspection of the DF of Hernquist models 
with Cuddeford anisotropy. 

(vi) To investigate the combined effect of the outer radial and in- 
ner tangential anisotropy, we performed an asymptotic analysis of 
one-component 71-7 models with Cuddeford anisotropy. We found 
that a steepening of the external density slope allows larger values 
of the central tangential anisotropy, independently of the value of 
the central density slope 7, thus confirming the hypothesis of a dy- 
namical interplay between the two regions of the models, and sup- 
porting the interpretation that Baes & Dejonghe (2002) proposed 
for Hernquist models. 

We notice that one of the major results of this study seems to 
be the generality of the density slope-anisotropy relation y(r) > 
2/3(r). It is natural to ask whether such density slope-anisotropy 
relation is even more general, i.e. it is an inequality necessarily 
obeyed by generic spherically symmetric, two-integrals systems 
with positive DF. At this stage we do not have a proof of this con- 
jecture, but we are not aware of any counter-example. Actually, 
we have additional evidences supporting this conjecture: for exam- 
ple Michele Trenti kindly provided us with a large set of numeri- 
cally computed f v models (Bertin & Trenti 2003), and all of them, 
without exception, satisfy the inequality f(r) > 2d(r) at all radii. 
Moreover, it is trivial to show that spherical systems in which the 
density can be written as p = A(r)/(^I'), with / monotonically in- 



creasing function of ty, all obey to y(r) > 2/3(r) when supported 
by a positive DF (see also the comment after equation IA9I ). We 
stress that these models do not belong to the family of generalized 
Cuddeford systems. Other distributions of orbital anisotropy that 
are not of the Cuddeford family (even though they could be approx- 
imated by specific choices of generalized Cuddeford distributions) 
have been reported by Mamon & Lokas (2005), Wojtak et al. (2008, 
who went further to show that also the DF was not of the OM or 
Cuddeford forms), Ascasibar et al. (2008) from the analysis of ha- 
los in cosmological simulations, or proposed in terms of specific 
DFs (e.g., Gerhard 1991; Louis 1993; Cuddeford & Louis 1995): 
it would be interesting to check the y(r) > 2B(r) inequality in 
these systems. In any case, we note that numerical simulations are 
known to produce correlations between (3 and 7 (e.g., see Hansen & 
Moore 2006, Mamon et al. 2006). We finally conclude by noticing 
that, if the inequality 7(7-) > 2/3(r) is universal (for spherical sys- 
tems), then the so-called mass-anisotropy degeneracy could be less 
severe than what is sometimes feared, as orbital anisotropy would 
be in some sense controlled by the local density slope of the stellar 
distribution in galaxies (in the inner regions where 7 < 2). This 
could be an important constraint in observational works. 
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APPENDIX A: THE DENSITY SLOPE-ANISOTROPY 
INEQUALITY FOR GENERALIZED CUDDEFORD 
SYSTEMS 

The radial and tangential velocity dispersion profiles of each den- 
sity component of a multi-component generalized Cuddeford sys- 
tem are given by 



2 / r 2 ,3 
pa r — I jv T a v 



47rB(r,a) / [2(* T - Q)] a+3/2 h(Q)dQ, (Al) 



2 / r 2 ,3 

Po- t = / fv t d v 



= 47rC*(r,a) / [2(* T - Q)] a+ ' A/2 h(Q)dQ, (A2) 
Jo 



where 



B(r,a) 



C(r,a) 



r(a + l) 
4 F(a + 5/2) (1 + r 2/ r 2.)a+ 



Err 



ypH I> + 2) 



2 r(a + 5/2) 51 (l + r 2 /^i)" +2 ' 

i 

Then, from equation $4$ one has 

P{ ' 2B(r,a)' 

and simple algebra proves equation d23t . 



WiT 



(A3) 
(A4) 

(A5) 
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We now show that the inequality ~f(r) > 2/3(r) holds at all 
radii in each density component of consistent generalized Cudde- 
ford systems. First, we relate the logarithmic density slope 7(7") = 
-dlnp/dln r to the NCi i|29} as follows: 



p dln(p/A) 
din r 



> ^£ = —£. = 
~ dr dr A rA 

so that the NCi can be simply rewritten as 

, . din A 
7(7-) > 



dlnr 



(A6) 



(A7) 



In other words, all consistent generalized Cuddeford systems sat- 
isfy equation JA7b at each radius. Now it is easy to verify that the 
functions A, B, and C, given in equations d24t . jA3\ , and (A4), 
satisfy the identity 



din .4 
din r 



= 2 



1 - 



C(r,a) 
2B{r,a) 



= 2/3 



(A8) 



for arbitrary a, Wi and r a i, so that equations dA7t and d!A8b show 
that the inequality 7(f) > 2/3(r) is just another way to express the 
NCi. Identity dA8t can be proved by elementary algebra: 

din A 



= -2a 
dlnr 

+2(a + l)^"' r2+ ^ 



(A9) 



2 W+2 



(1+r 2 /r 2 i)a4 



where in each term of the sum we added and subtracted u>ir ai . Sim- 
plification and comparison with equation l |23t conclude the proof. 

We note that identity jAS\ is actually a special case of a more 
general result reported in Baes & Dejonghe (2002) and Baes & van 
Hese (2007), which holds for all spherical systems whose DF, after 
integration over velocity space, leads to the factorization p(r) = 
A(r)f(ty); for such systems (that also include our generalized 
Cuddeford systems) it can be proved that 2/3(r) — — din A/d\n r. 
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